A model of alcoholic liver disease based on different hepatotoxics leading to liver cancer.

The worst-case scenario related to alcoholic liver disease (ALD) arises after a long period of exposure to the harmful effect of alcohol consumption along with other hepatotoxics. ALD encompasses a broad spectrum of liver-associated disorders, such as steatosis, steatohepatitis, fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). Based on the chronic administration of different hepatotoxics, including ethanol, sucrose, lipopolysaccharide, and low doses of diethylnitrosamine over a short period, here we aimed to develop a multiple hepatotoxic (MHT)-ALD model in the mouse that recapitulates the human ALD-associated disorders. We demonstrated that the MHT-ALD model induces ADH1A and NXN, an ethanol metabolizer and a redox-sensor enzyme, respectively; promotes steatosis associated with the induction of the lipid droplet forming FSP27, inflammation identified by the infiltration of hepatic neutrophils-positive to LY-6G marker, and the increase of MYD88 level, a protein involved in inflammatory response; and stimulates the early appearance of cellular senescence identified by the senescence markers SA-ß-gal activity and p-H2A.XSer139. It also induces fibrosis associated with increased desmin, a marker of hepatic stellate cells whose activation leads to the deposition of collagen fibers, accompanied by cell death and compensatory proliferation revealed by increased CASP3-mediated apoptosis, and KI67- and PCNA-proliferation markers, respectively. It also induces histopathological traits of malignancy and the level of the HCC marker, GSTP1. In conclusion, we provide a useful model for exploring the chronological ALD-associated alterations and stages, and addressing therapeutic approaches.

Mobility networks in Greater Mexico City.

Based on more than 11 billion geolocated cell phone records from 33 million different devices, daily mobility networks were constructed over a 15-month period for Greater Mexico City, one of the largest and most diverse metropolitan areas globally. The time frame considered spans the entire year of 2020 and the first three months of 2021, enabling the analysis of population movement dynamics before, during, and after the COVID-19 health contingency. The nodes within the 456 networks represent the basic statistical geographic areas (AGEBs) established by the National Institute of Statistics, Geography, and Informatics (INEGI) in Mexico. This framework facilitates the integration of mobility data with numerous indicators provided by INEGI. Edges connecting these nodes represent movement between AGEBs, with edge weights indicating the volume of trips from one AGEB to another. This extensive dataset allows researchers to uncover travel patterns, cross-reference data with socio-economic indicators, and conduct segregation studies, among other potential analyses.

Computational single cell oncology: state of the art.

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

Exploring disparities and novel insights into metabolic and cardiovascular comorbidities among COVID-19 patients in Mexico.

Introduction: The COVID-19 pandemic, especially its early stages, sparked extensive discussions regarding the potential impact of metabolic and cardiovascular comorbidities on the severity and fatality of SARS-CoV-2 infection, yielding inconclusive outcomes. In this study, we delve into the prevalence of metabolic and cardiovascular comorbidities within COVID-19 patients in Mexico. Methods: Employing a retrospective observational study design, we collected data from official databases encompassing COVID-19 patients admitted to both public and private hospitals in Mexico City. Results: Our investigation unveiled a noteworthy incongruity in the prevalence of metabolic and cardiovascular comorbidities among COVID-19 patients, with a particular emphasis on obesity, hypertension, and diabetes. This incongruity manifests as location-dependent phenomena, where the prevalence of these comorbidities among COVID-19 patients significantly deviates from the reported values for the general population in each specific location. Discussion: These findings underscore the critical importance of screening for metabolic and cardiovascular comorbidities in COVID-19 patients and advocate for the necessity of tailored interventions for this specific population. Furthermore, our study offers insights into the intricate interplay between COVID-19 and metabolic and cardiovascular comorbidities, serving as a valuable foundation for future research endeavors and informing clinical practice.

Differential effects of follicle-stimulating hormone glycoforms on the transcriptome profile of cultured rat granulosa cells as disclosed by RNA-seq.

It has been documented that variations in glycosylation on glycoprotein hormones, confer distinctly different biological features to the corresponding glycoforms when multiple in vitro biochemical readings are analyzed. We here applied next generation RNA sequencing to explore changes in the transcriptome of rat granulosa cells exposed for 0, 6, and 12 h to 100 ng/ml of four highly purified follicle-stimulating hormone (FSH) glycoforms, each exhibiting different glycosylation patterns: human pituitary FSH18/21 and equine FSH (eqFSH) (hypo-glycosylated), and human FSH24 and chinese-hamster ovary cell-derived human recombinant FSH (recFSH) (fully-glycosylated). Total RNA from triplicate incubations was prepared from FSH glycoform-exposed cultured granulosa cells obtained from DES-pretreated immature female rats, and RNA libraries were sequenced in a HighSeq 2500 sequencer (2 × 125 bp paired-end format, 10-15 × 106 reads/sample). The computational workflow focused on investigating differences among the four FSH glycoforms at three levels: gene expression, enriched biological processes, and perturbed pathways. Among the top 200 differentially expressed genes, only 4 (0.6%) were shared by all 4 glycoforms at 6 h, whereas 118 genes (40%) were shared at 12 h. Follicle-stimulating hormone glycocoforms stimulated different patterns of exclusive and associated up regulated biological processes in a glycoform and time-dependent fashion with more shared biological processes after 12 h of exposure and fewer treatment-specific ones, except for recFSH, which exhibited stronger responses with more specifically associated processes at this time. Similar results were found for down-regulated processes, with a greater number of processes at 6 h or 12 h, depending on the particular glycoform. In general, there were fewer downregulated than upregulated processes at both 6 h and 12 h, with FSH18/21 exhibiting the largest number of down-regulated associated processes at 6 h while eqFSH exhibited the greatest number at 12 h. Signaling cascades, largely linked to cAMP-PKA, MAPK, and PI3/AKT pathways were detected as differentially activated by the glycoforms, with each glycoform exhibiting its own molecular signature. These data extend previous observations demonstrating glycosylation-dependent differential regulation of gene expression and intracellular signaling pathways triggered by FSH in granulosa cells. The results also suggest the importance of individual FSH glycoform glycosylation for the conformation of the ligand-receptor complex and induced signalling pathways.

The Role of Carbohydrate Intake on the Gut Microbiome: A Weight of Evidence Systematic Review.

(1) Background: Carbohydrates are the most important source of nutritional energy for the human body. Carbohydrate digestion, metabolism, and their role in the gut microbiota modulation are the focus of multiple studies. The objective of this weight of evidence systematic review is to investigate the potential relationship between ingested carbohydrates and the gut microbiota composition at different taxonomic levels. (2) Methods: Weight of evidence and information value techniques were used to evaluate the relationship between dietary carbohydrates and the relative abundance of different bacterial taxa in the gut microbiota. (3) Results: The obtained results show that the types of carbohydrates that have a high information value are: soluble fiber with Bacteroides increase, insoluble fiber with Bacteroides and Actinobacteria increase, and Firmicutes decrease. Oligosaccharides with Lactobacillus increase and Enterococcus decrease. Gelatinized starches with Prevotella increase. Starches and resistant starches with Blautia decrease and Firmicutes increase. (4) Conclusions: This work provides, for the first time, an integrative review of the subject by using statistical techniques that have not been previously employed in microbiota reviews.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico.

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Drug repurposing for Basal breast cancer subpopulations using modular network signatures.

Breast cancer is characterized as being a heterogeneous pathology with a broad phenotype variability. Breast cancer subtypes have been developed in order to capture some of this heterogeneity. Each of these breast cancer subtypes, in turns retains varied characteristic features impacting diagnostic, prognostic and therapeutics. Basal breast tumors, in particular have been challenging in these regards. Basal breast cancer is often more aggressive, of rapid evolution and no tailor-made targeted therapies are available yet to treat it. Arguably, epigenetic variability is behind some of these intricacies. It is possible to further classify basal breast tumor in groups based on their non-coding transcriptome and methylome profiles. It is expected that these groups will have differences in survival as well as in sensitivity to certain classes of drugs. With this in mind, we implemented a computational learning approach to infer different subpopulations of basal breast cancer (from TCGA multi-omic data) based on their epigenetic signatures. Such epigenomic signatures were associated with different survival profiles; we then identified their associated gene co-expression network structure, extracted a signature based on modules within these networks, and use these signatures to find and prioritize drugs (in the LINCS dataset) that may be used to target these types of cancer. In this way we are introducing the analytical workflow for an epigenomic signature-based drug repurposing structure.

COVID-19 reinfections in Mexico City: implications for public health.

Background: Since its appearance, COVID-19 has immensely impacted our society. Public health measures, from the initial lockdowns to vaccination campaigns, have mitigated the crisis. However, SARS-CoV-2's persistence and evolving variants continue to pose global threats, increasing the risk of reinfections. Despite vaccination progress, understanding reinfections remains crucial for informed public health responses. Methods: We collected available data on clinical and genomic information for SARS-CoV-2 samples from patients treated in Mexico City from 2020 epidemiological week 10 to 2023 epidemiological week 06 encompassing the whole public health emergency's period. To identify clinical data we utilized the SISVER (Respiratory Disease Epidemiological Surveillance System) database for SARS-CoV-2 patients who received medical attention in Mexico City. For genomic surveillance we analyzed genomic data previously uploaded to GISAID generated by Mexican institutions. We used these data sources to generate descriptors of case number, hospitalization, death and reinfection rates, and viral variant prevalence throughout the pandemic period. Findings: The fraction of reinfected individuals in the COVID-19 infected population steadily increased as the pandemic progressed in Mexico City. Most reinfections occurred during the fifth wave (40%). This wave was characterized by the coexistence of multiple variants exceeding 80% prevalence; whereas all other waves showed a unique characteristic dominant variant (prevalence >95%). Shifts in symptom patient care type and severity were observed, 2.53% transitioned from hospitalized to ambulatory care type during reinfection and 0.597% showed the opposite behavior; also 7.23% showed a reduction in severity of symptoms and 6.05% displayed an increase in severity. Unvaccinated individuals accounted for the highest percentage of reinfections (41.6%), followed by vaccinated individuals (31.9%). Most reinfections occurred after the fourth wave, dominated by the Omicron variant; and after the vaccination campaign was already underway. Interpretation: Our analysis suggests reduced infection severity in reinfections, evident through shifts in symptom severity and care patterns. Unvaccinated individuals accounted for most reinfections. While our study centers on Mexico City, its findings may hold implications for broader regions, contributing insights into reinfection dynamics.

Probability of hospitalisation and death among COVID-19 patients with comorbidity during outbreaks occurring in Mexico City.

Background: We compared the probability of hospitalization and death caused by COVID-19 in patients with comorbidities during three periods defined for this study: first-wave (FW), interwave period (IP), and second-wave (SW) observed in Mexico City. Methods: In this registry-based study, we included individuals over 20 years of age. During the FW (symptomatic), the IP, and the SW (symptomatic and asymptomatic), participants were diagnosed using nasopharyngeal swabs. Symptomatic individuals with risk factors for serious disease or death were referred to the hospital. SARS-CoV-2 infection was defined by RT-qPCR in all hospitalized patients. All data were added to the SISVER database. Bayesian analysis and False Discovery Rate were used for further evaluation. Results: The study included 2 260 156 persons (mean age of 43.1 years). Of these, 8.6% suffered from DM, 11.6% arterial hypertension, and 9.7% obesity. Of the total, 666 694 persons tested positive (29.5%). Of the infected persons, a total of 85 587 (12.8%) were hospitalized: 24 023 in the FW; 16 935 in the IP, and 44 629 in the SW. Of the hospitalized patients, there were 42 979 deaths (50.2%), in the FW, 11 964 (49.8%), in the IP, 6794 (40.1%), and in the SW 24 221 (54.3%). The probability of death among individuals hospitalized with or without comorbidities increased consistently in all age groups. A significant increase in the Fatality Rate was observed in individuals with comorbidities (1.36E-19< = FDR< = 3.36E-2). A similar trend was also observed in individuals without comorbidities (1.03E-44< = FDR< = 5.58E-4). Conclusions: The data from this study show a considerable increase in the number of detected cases of infection between the FW and SW. In addition, 12.8% of those infected were hospitalized for severe COVID-19. A high mortality rate was observed among hospitalized patients (>50%). An age-dependent probability of death was observed with a positive trend in hospitalized patients with and without comorbidities.

Modular reactivation of Mexico City after COVID-19 lockdown.

BACKGROUND: During the COVID-19 pandemic, the slope of the epidemic curve in Mexico City has been quite unstable. Changes in human activity led to changes in epidemic activity, hampering attempts at economic and general reactivation of the city. METHODS: We have predicted that where a fraction of the population above a certain threshold returns to the public space, the negative tendency of the epidemic curve will revert. Such predictions were based on modeling the reactivation of economic activity after lockdown using an epidemiological model resting upon a contact network of Mexico City derived from mobile device co-localization. We modeled scenarios with different proportions of the population returning to normalcy. Null models were built using the Jornada Nacional de Sana Distancia (the Mexican model of elective lockdown). There was a mobility reduction of 75% and no mandatory mobility restrictions. RESULTS: We found that a new peak of cases in the epidemic curve was very likely for scenarios in which more than 5% of the population rejoined the public space. The return of more than 50% of the population synchronously will unleash a magnitude similar to the one predicted with no mitigation strategies. By evaluating the tendencies of the epidemic dynamics, the number of new cases registered, hospitalizations, and recent deaths, we consider that reactivation following only elective measures may not be optimal under this scenario. CONCLUSIONS: Given the need to resume economic activities, we suggest alternative measures that minimize unnecessary contacts among people returning to the public space. We evaluated that "encapsulating" reactivated workers (that is, using measures to reduce the number of contacts beyond their influential community in the contact network) may allow reactivation of a more significant fraction of the population without compromising the desired tendency in the epidemic curve.

Introduction to Genomic Network Reconstruction for Cancer Research.

High-throughput genomic technologies have revolutionized the study of cancer. Current research in oncology is now limited more for the capacity of analyzing and interpreting data, rather than the availability of data itself. Integrative approaches to obtain functional information from data are at the core of the disciplines gathered under the systems biology banner. In this context, network models have been used to study cancer, from the identification of key molecules involved in the disease to the discovery of functional alterations associated with specific manifestations of the disease.In this chapter, we describe the state of the art of network reconstruction from genomic data, with an emphasis in gene expression experiments. We explore the strengths and limitations of correlation, Bayesian, and information theoretic approaches to network reconstruction. We present tools that leverage the flexibility of network science to gain a deeper understanding of cancer biology.

Early Genomic, Epidemiological, and Clinical Description of the SARS-CoV-2 Omicron Variant in Mexico City.

Omicron is the most mutated SARS-CoV-2 variant-a factor that can affect transmissibility, disease severity, and immune evasiveness. Its genomic surveillance is important in cities with millions of inhabitants and an economic center, such as Mexico City. Results. From 16 November to 31 December 2021, we observed an increase of 88% in Omicron prevalence in Mexico City. We explored the R346K substitution, prevalent in 42% of Omicron variants, known to be associated with immune escape by monoclonal antibodies. In a phylogenetic analysis, we found several independent exchanges between Mexico and the world, and there was an event followed by local transmission that gave rise to most of the Omicron diversity in Mexico City. A haplotype analysis revealed that there was no association between haplotype and vaccination status. Among the 66% of patients who have been vaccinated, no reported comorbidities were associated with Omicron; the presence of odynophagia and the absence of dysgeusia were significant predictor symptoms for Omicron, and the RT-qPCR Ct values were lower for Omicron. Conclusions. Genomic surveillance is key to detecting the emergence and spread of SARS-CoV-2 variants in a timely manner, even weeks before the onset of an infection wave, and can inform public health decisions and detect the spread of any mutation that may affect therapeutic efficacy.

Territorial Strategy of Medical Units for Addressing the First Wave of the COVID-19 Pandemic in the Metropolitan Area of Mexico City: Analysis of Mobility, Accessibility and Marginalization.

BACKGROUND: The COVID-19 pandemic has caused an exponential increase in the demand for medical care worldwide. In Mexico, the COVID Medical Units (CMUs) conversion strategy was implemented. OBJECTIVE: To evaluate the CMU coverage strategy in the Mexico City Metropolitan Area (MCMA) by territory. MATERIALS: The CMU directory was used, as were COVID-19 infection and mobility statistics and Mexican 2020 census information at the urban geographic area scale. The degree of urban marginalization by geographic area was also considered. METHOD: Using descriptive statistics and the calculation of a CMU accessibility index, population aggregates were counted based on coverage radii. In addition, two regression models are proposed to explain (1) the territorial and temporal trend of COVID-19 infections in the MCMA and (2) the mobility of the COVID-infected population visiting medical units. RESULTS: The findings of the evaluation of the CMU strategy were (1) in the MCMA, COVID-19 followed a pattern of contagion from the urban center to the periphery; (2) given the growth in the number of cases and the overload of medical units, the population traveled greater distances to seek medical care; (3) after the CMU strategy was evaluated at the territory level, it was found that 9 out of 10 inhabitants had a CMU located approximately 7 km away; and (4) at the metropolitan level, the lowest level of accessibility to the CMU was recorded for the population with the highest levels of marginalization, i.e., those residing in the urban periphery.

The Evolutionary Landscape of SARS-CoV-2 Variant B.1.1.519 and Its Clinical Impact in Mexico City.

The SARS-CoV-2 pandemic is one of the most concerning health problems around the globe. We reported the emergence of SARS-CoV-2 variant B.1.1.519 in Mexico City. We reported the effective reproduction number (Rt) of B.1.1.519 and presented evidence of its geographical origin based on phylogenetic analysis. We also studied its evolution via haplotype analysis and identified the most recurrent haplotypes. Finally, we studied the clinical impact of B.1.1.519. The B.1.1.519 variant was predominant between November 2020 and May 2021, reaching 90% of all cases sequenced in February 2021. It is characterized by three amino acid changes in the spike protein: T478K, P681H, and T732A. Its Rt varies between 0.5 and 2.9. Its geographical origin remain to be investigated. Patients infected with variant B.1.1.519 showed a highly significant adjusted odds ratio (aOR) increase of 1.85 over non-B.1.1.519 patients for developing a severe/critical outcome (p = 0.000296, 1.33-2.6 95% CI) and a 2.35-fold increase for hospitalization (p = 0.005, 1.32-4.34 95% CI). The continuous monitoring of this and other variants will be required to control the ongoing pandemic as it evolves.

Automated Reverse Transcription Polymerase Chain Reaction Data Analysis for Sars-CoV-2 Detection.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic is a current public health concern. Rapid diagnosis is crucial, and reverse transcription polymerase chain reaction (RT-PCR) is presently the reference standard for SARS-CoV-2 detection. OBJECTIVE: Automated RT-PCR analysis (ARPA) is a software designed to analyze RT-PCR data for SARSCoV-2 detection. ARPA loads the RT-PCR data, classifies each sample by assessing its amplification curve behavior, evaluates the experiment's quality, and generates reports. METHODS: ARPA was implemented in the R language and deployed as a Shiny application. We evaluated the performance of ARPA in 140 samples. The samples were manually classified and automatically analyzed using ARPA. RESULTS: ARPA had a true-positive rate = 1, true-negative rate = 0.98, positive-predictive value = 0.95, and negative-predictive value = 1, with 36 samples correctly classified as positive, 100 samples correctly classified as negative, and two samples classified as positive even when labeled as negative by manual inspection. Two samples were labeled as invalid by ARPA and were not considered in the performance metrics calculation. CONCLUSIONS: ARPA is a sensitive and specific software that facilitates the analysis of RT-PCR data, and its implementation can reduce the time required in the diagnostic pipeline.

Highly connected, non-redundant microRNA functional control in breast cancer molecular subtypes.

Breast cancer is a complex, heterogeneous disease at the phenotypic and molecular level. In particular, the transcriptional regulatory programs are known to be significantly affected and such transcriptional alterations are able to capture some of the heterogeneity of the disease, leading to the emergence of breast cancer molecular subtypes. Recently, it has been found that network biology approaches to decipher such abnormal gene regulation programs, for instance by means of gene co-expression networks, have been able to recapitulate the differences between breast cancer subtypes providing elements to further understand their functional origins and consequences. Network biology approaches may be extended to include other co-expression patterns, like those found between genes and non-coding transcripts such as microRNAs (miRs). As is known, miRs play relevant roles in the establishment of normal and anomalous transcription processes. Commodore miRs (cdre-miRs) have been defined as miRs that, based on their connectivity and redundancy in co-expression networks, are potential control elements of biological functions. In this work, we reconstructed miR-gene co-expression networks for each breast cancer molecular subtype, from high throughput data in 424 samples from the Cancer Genome Atlas consortium. We identified cdre-miRs in three out of four molecular subtypes. We found that in each subtype, each cdre-miR was linked to a different set of associated genes, as well as a different set of associated biological functions. We used a systematic literature validation strategy, and identified that the associated biological functions to these cdre-miRs are hallmarks of cancer such as angiogenesis, cell adhesion, cell cycle and regulation of apoptosis. The relevance of such cdre-miRs as actionable molecular targets in breast cancer is still to be determined from functional studies.

In silico identification of excretory/secretory proteins and drug targets in monogenean parasites.

The Excretory/Secretory (ES) proteins of parasites are involved in invasion and colonization of their hosts. In addition, since ES proteins circulate in the extracellular space, they can be more accessible to drugs than other proteins, which makes ES proteins optimal targets for the development of new and better pharmacological strategies. Monogeneans are a group of parasitic Platyhelminthes that includes some pathogenic species problematic for finfish aquaculture. In the present study, 8297 putative ES proteins from four monogenean species which genomic resources are publicly available were identified and functionally annotated by bioinformatic tools. Additionally, for comparative purposes, ES proteins in other parasitic and free-living platyhelminths were identified. Based on data from the monogenean Gyrodactylus salaris, 15 ES proteins are considered potential drug targets. One of them showed homology to 10 cathepsins with known 3D structure. A docking molecular analysis uncovered that the anthelmintic emodepside shows good affinity to these cathepsins suggesting that emodepside can be experimentally tested as a monogenean's cathepsin inhibitor.

An Information Theoretical Multilayer Network Approach to Breast Cancer Transcriptional Regulation.

Breast cancer is a complex, highly heterogeneous disease at multiple levels ranging from its genetic origins and molecular processes to clinical manifestations. This heterogeneity has given rise to the so-called intrinsic or molecular breast cancer subtypes. Aside from classification, these subtypes have set a basis for differential prognosis and treatment. Multiple regulatory mechanisms-involving a variety of biomolecular entities-suffer from alterations leading to the diseased phenotypes. Information theoretical approaches have been found to be useful in the description of these complex regulatory programs. In this work, we identified the interactions occurring between three main mechanisms of regulation of the gene expression program: transcription factor regulation, regulation via noncoding RNA, and epigenetic regulation through DNA methylation. Using data from The Cancer Genome Atlas, we inferred probabilistic multilayer networks, identifying key regulatory circuits able to (partially) explain the alterations that lead from a healthy phenotype to different manifestations of breast cancer, as captured by its molecular subtype classification. We also found some general trends in the topology of the multi-omic regulatory networks: Tumor subtype networks present longer shortest paths than their normal tissue counterpart; epigenomic regulation has frequently focused on genes enriched for certain biological processes; CpG methylation and miRNA interactions are often part of a regulatory core of conserved interactions. The use of probabilistic measures to infer information regarding theoretical-derived multilayer networks based on multi-omic high-throughput data is hence presented as a useful methodological approach to capture some of the molecular heterogeneity behind regulatory phenomena in breast cancer, and potentially other diseases.

Gene-Microbiome Co-expression Networks in Colon Cancer.

It is known that cancer onset and development arise from complex, multi-factorial phenomena spanning from the molecular, functional, micro-environmental, and cellular up to the tissular and organismal levels. Important advances have been made in the systematic analysis of the molecular (mostly genomic and transcriptomic) within large studies of high throughput data such as The Cancer Genome Atlas collaboration. However, the role of the microbiome in the induction of biological changes needed to reach these pathological states remains to be explored, largely because of scarce experimental data. In recent work a non-standard bioinformatics strategy was used to indirectly quantify microbial abundance from TCGA RNA-seq data, allowing the evaluation of the microbiome in well-characterized cancer patients, thus opening the way to studies incorporating the molecular and microbiome dimensions altogether. In this work, we used such recently described approaches for the quantification of microbial species alongside with gene expression. With this, we will reconstruct bipartite networks linking microbial abundance and gene expression in the context of colon cancer, by resorting to network reconstruction based on measures from information theory. The rationale is that microbial communities may induce biological changes important for the cancerous state. We analyzed changes in microbiome-gene interactions in the context of early (stages I and II) and late (stages III and IV) colon cancer, studied changes in network descriptors, and identify key discriminating features for early and late stage colon cancer. We found that early stage bipartite network is associated with the establishment of structural features in the tumor cells, whereas late stage is related to more advance signaling and metabolic features. This functional divergence thus arise as a consequence of changes in the organization of the corresponding gene-microorganism co-expression networks.